home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlasq2.z / dlasq2
Encoding:
Text File  |  2002-10-03  |  4.3 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSQQQQ2222((((3333SSSS))))                                                          DDDDLLLLAAAASSSSQQQQ2222((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASQ2 - compute all the eigenvalues of the symmetric positive definite
  10.      tridiagonal matrix associated with the qd array Z to high relative
  11.      accuracy are computed to high relative accuracy, in the absence of
  12.      denormalization, underflow and overflow
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  15.      SUBROUTINE DLASQ2( N, Z, INFO )
  16.  
  17.          INTEGER        INFO, N
  18.  
  19.          DOUBLE         PRECISION Z( * )
  20.  
  21. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  22.      These routines are part of the SCSL Scientific Library and can be loaded
  23.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  24.      directs the linker to use the multi-processor version of the library.
  25.  
  26.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  27.      4 bytes (32 bits). Another version of SCSL is available in which integers
  28.      are 8 bytes (64 bits).  This version allows the user access to larger
  29.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  30.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  31.      only one of the two versions; 4-byte integer and 8-byte integer library
  32.      calls cannot be mixed.
  33.  
  34. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  35.      DLASQ2 computes all the eigenvalues of the symmetric positive definite
  36.      tridiagonal matrix associated with the qd array Z to high relative
  37.      accuracy are computed to high relative accuracy, in the absence of
  38.      denormalization, underflow and overflow. To see the relation of Z to the
  39.      tridiagonal matrix, let L be a unit lower bidiagonal matrix with
  40.      subdiagonals Z(2,4,6,,..) and let U be an upper bidiagonal matrix with
  41.      1's above and diagonal Z(1,3,5,,..). The tridiagonal is L*U or, if you
  42.      prefer, the symmetric tridiagonal to which it is similar.
  43.  
  44.      Note : DLASQ2 defines a logical variable, IEEE, which is true on machines
  45.      which follow ieee-754 floating-point standard in their handling of
  46.      infinities and NaNs, and false otherwise. This variable is passed to
  47.      DLASQ3.
  48.  
  49.  
  50. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  51.      N     (input) INTEGER
  52.            The number of rows and columns in the matrix. N >= 0.
  53.  
  54.      Z     (workspace) DOUBLE PRECISION array, dimension ( 4*N )
  55.            On entry Z holds the qd array. On exit, entries 1 to N hold the
  56.            eigenvalues in decreasing order, Z( 2*N+1 ) holds the trace, and Z(
  57.            2*N+2 ) holds the sum of the eigenvalues. If N > 2, then Z( 2*N+3 )
  58.            holds the iteration count, Z( 2*N+4 ) holds NDIVS/NIN^2, and Z(
  59.            2*N+5 ) holds the percentage of shifts that failed.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAASSSSQQQQ2222((((3333SSSS))))                                                          DDDDLLLLAAAASSSSQQQQ2222((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      INFO  (output) INTEGER
  75.            = 0: successful exit
  76.            < 0: if the i-th argument is a scalar and had an illegal value,
  77.            then INFO = -i, if the i-th argument is an array and the j-entry
  78.            had an illegal value, then INFO = -(i*100+j) > 0: the algorithm
  79.            failed = 1, a split was marked by a positive value in E = 2,
  80.            current block of Z not diagonalized after 30*N iterations (in inner
  81.            while loop) = 3, termination criterion of outer while loop not met
  82.            (program created more than N unreduced blocks)
  83.  
  84. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  85.      The shifts are accumulated in SIGMA. Iteration count is in ITER.  Ping-
  86.      pong is controlled by PP (alternates between 0 and 1).
  87.  
  88.  
  89. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  90.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  91.  
  92.      This man page is available only online.
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.